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A numerical study of bifurcations 
in a barotropic shear flow 
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1. Introduction 

In the last few years, more and more evidence has emerged suggesting that tran- 
sition to turbulence may be viewed as a succession of bifurcations to deterministic 
chaos. Most experimental and numerical observations have, for the most part, been 
restricted to Rayleigh-Benard convection and Taylor-Couette flow between concen- 
tric cylinders. The goal of the present study is to accurately describe the bifurcation 
sequence leading to chaos in a two-dimensional temporal free shear layer on the 0- 
plane. 

The /?-plane is a locally- Cartesian reduction of the equations describing the dy- 
namics of a shallow layer of fluid on a rotating spherical planet. It is a valid model 
for large-scale flows of interest in meteorology and oceanography. The mathematical 
formulation involves the following vorticity equation 

(Ji + V 2 * + «7(V 2 ¥,¥) + (/? - 

where U(y) is the basic hyperbolic-tangent velocity profile, 9 is the perturbation 
stream function, and Re is the Reynolds number. The symbol 0 denotes the gradient 
of the planetary vorticity. Inviscid linearized stability analyses (Dickinson and Clare 
1973) reveal that the flow is unstable within a neutral curve in the 0 — k plane as 
sketched in Figure 1 (from Burns & Maslowe 1983), k denoting the streamwise 
wavenumber. Thus 0 is a convenient control parameter with which one can bring 
about qualitative changes in the attractor as a function of the supercriticality 
\0 — 0 C \' The limit of the classical mixing layer is approached as 0 goes to zero. 

2. Numerical approach 

The three-dimensional code developed by M. Rogers has been adapted to the 
present problem by adding the appropriate 0 terms in the equations. The basic 
velocity profile has been changed to U(y) = tanh(t/) (a body force is applied to hold 
that velocity profile against viscous diffusion). We recall that periodic boundary 
conditions are chosen in the streamwise direction and that a Cain mapping is used in 
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the cross-stream direction. The exponential tails in y for the perturbation variables 
are known to decay more slowly than in the (3 = 0 case and the value of the B 
stretching factor in y had to be increased up to 10 (typical values are around 2) to 
properly account for this phenomenon. The numerical code is well-behaved with 16 
by 96 grids, but the execution is slow (8 Mflops on the Cray 2) due to the small 
length of the vectors. 

3. The first Hopf bifurcation 

Numerical simulations were conducted at Re = 10 00 a nd a streamwise box length 
corresponding to the critical wavenumber k c = y/2/3. The critical value f3 c — 
4/(3\/3) was well approximated and linear growth rates were also found to be 
consistent with inviscid stability theory. 

As j3 decreases below /? c , the basic flow bifurcates to a finite-amplitude traveling - 
wave state which takes the form of a vortex as shown on Figure 2. As a result of the 
/? effect, the Kelvin’s cat’s-eye pattern is shifted off the axis y = 0. The traveling 
wave is characterized by a circular frequency <j) and normalized amplitude 

|A| = max v(y,k c ). 
y 

The normal form pertaining to a supercritical Hopf bifurcation is known to be 

^ = -i[u e + «„(/* - 13c)] A - T\A\ 2 A, (1) 

where A is the complex amplitude of the wave and u; c the frequency at k = k c , 
(3 = /? c . Other quantities appearing in (1) are = d(v/dfi\ c and the Landau 
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FIGURE 2. Isovorticity contours at 0 = 0.7, Re = 1000, k = k c . 

constant T = T r + il\. An elementary analysis of (1) indicates that as 0 decreases 
below 0 cy the solution of (1) bifurcates from the basic flow A = 0 to a limit cycle 
of amplitude 

\A\ 2 oc 0c — (3 (2) 

and frequency 

v - U c oc Pc - 0- (3) 

The quantities \A\ 2 and u> are represented as a function of 0 on Figures 3 and 4 
respectively. The exact linear value of u) has also been plotted on Figure 4. One 
notes that \A\ 2 and a? do vary linearly with decreasing 0 in the vicinity of 0 C) in 
agreement with (2) and (3). The value of u; at 0 = 0 C is indeed equal to u) c = y/2/3. 
Finally one notes a nonlinear frequency shift in the value of w when compared to 
linear estimates. 

We attempted to compare the values of the slopes of these curves with the analyt- 
ical results obtained by Churilov (1988), Churilov and Shukhman (1986, 1987), and 
Lipps (1965). There were large discrepancies and we suspect that the critical-layer 
structure is not satisfactorily resolved with the present distribution of points in the 
y direction. 

4. A preliminary look at pairing in the presence of dispersion 

To test the validity of the code, numerical simulations of the evolution of 3 vor- 
tices were conducted in a computational domain of length 6n/k c at 0 = 0.2. The 



Figure 4. u versus /3 at Re = 1000, k = k c . — , nonlinear frequency; , linear 

frequency. 
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fundamental Jb c was given a finite energy level at t = 0, the 2 subharmonics at 
k c /3 and 2Jb c /3 being initialized to low levels. Figures 5 a, b, c show the resulting 
evolution. One notes the nutation of the vortical structures in the course of time, 
indicated by a tilting motion. Pairing eventually takes place between 490 and 500 
time steps with a corresponding energy transfer from the k c to the 2k c /3 compo- 
nent. At (3 = 0.2, the k c /3 component lies in the linearly stable domain and it is not 
triggered by nonlinear interactions. Thus, no further pairing events are observed 
and the asymptotic state consists of 2 vortices. This is distinctly different from the 
homogeneous mixing layer case where pairings take place until one observes a triple 
vortex condition in the computational domain. 

5. Future work 

The results obtained so far, namely the fact that a Hopf bifurcation takes place 
at /? = /3 cy suggest that nontrivial chaotic dynamics are likely to occur within 
subregions of the neutral curve. To back up this claim, we note that in the vicinity 
of fle, the streamwise modulations of the complex amplitude A(x,f) are governed 
by the p.d.e. equivalent of (1), i.e.: 

dA . r 0 . dA i d 2 A . 2| 

— = - t [u e + Up (0 - &)] A - w k + -Wkk -^-T\A\A. 

This Ginzburg-Landau equation is known to give rise to chaos via a modulational 
instability mechanism [Keefe (1985)]. In the present context, we therefore suspect 
that amplitude and phase modulations applied to a row of vortices would induce 
a disordered evolution as the streamwise extent of the computational domain is 
increased. 

To test this hypothesis, we need to considerably increase the reliability and speed 
of the numerical code for low- wavenumber simulations. 
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